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Abstract. We prove here well-posedness and convergence to equilibria for the solu- 
tion trajectories associated to a model for solidification of a liquid content of a rigid 
container in a gravity field. We observe that the gravity effects, which can be neglected 
without considerable changes of the process on finite time intervals, have a substantial 
influence on the long time behavior of the evolution system. Without gravity, we find 
a temperature interval, in which all phase distributions with a prescribed total liquid 
contents are admissible equilibria, while, under the influence of gravity, the only equi- 
librium distribution in a connected container consists in two pure phases separated by 
one plane interface perpendicular to the gravity force. 
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1 Introduction 

In this paper we continue the discussion started in [7j, extending the model to the case 
in which also gravity effects are considered during the phase transition process. We 
derive a model for solid-liquid phase transition of a medium inside a rigid container. 
The main goal is to give a qualitative and quantitative description of the interaction 
between volume, pressure, phase, and temperature changes in the situation that the 
specific volume of the solid phase exceeds the specific volume of the liquid phase. We 
observe, in particular, that the solidification may take place at a temperature slightly 
above the critical temperature 9c ■ The overheating is due to the fact that the pressure 
decreases from the bottom to the top. A quantitative description of this phenomenon 
is given in here by means of the so-called Clausius-Clapeyron formula (cf. f l2.40p ). 

There is an abundant classical literature on the study of phase transition pro- 
cesses, see e.g. the monographs P, [3], [9] and the references therein. In particular, in 
[1] , the authors proposed to interpret a phase transition process in terms of a balance 
equation for macroscopic motions, and to include the possibility of voids, while the 
microscopic approach has been pursued in [5] in the case of two different densities gi 
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and ^2 for the two substances undergoing phase transitions. Let us, however, refer 
to the Introduction of [7j for a more detail description of the previous works in the 
hterature on this topic. 

The forces occurring as a result of solid-liquid phase transitions in small contain- 
ers are very strong, much stronger than gravity forces. In a bottle of water of less than 
one meter height for example, they differ by at least four orders of magnitude. From 
the quantitative viewpoint, the gravity effects can thus be neglected without consider- 
able changes of the process on finite time intervals. They have, however, a substantial 
influence on the long time behavior of the evolution system. Without gravity, we ob- 
serve a temperature interval, in which all phase distributions with a prescribed total 
liquid contents are admissible equilibria (cf. also [7j). If even a weak gravity field is 
assumed to be present, then the only equilibrium distribution in a connected con- 
tainer consists in two pure phases separated by one plane interface perpendicular to 
the gravity force. 

Here we proceed as follows: in Section [21 we derive a model describing the evo- 
lution of the process which is driven by an energy balance, a quasistatic momentum 
balance, and a phase dynamics equation. Still in Section [21 we verify the thermody- 
namic consistency of the model, and we study the equilibria. 

The well-posedness of the corresponding system of evolution equations is proved 
in Section [31 while and the study of the long-time behavior of solutions and convergence 
to equilibria is proved in the last Section [H 

2 The model 

As reference state, we consider a liquid substance contained in a bounded connected 
bottle f2 C with boundary of class C^'^ . The state variables are the absolute 
temperature 6 > 0, the displacement u G M^, and the phase variable x ^ [0, 1] • The 
value X = means solid, x = 1 means liquid, x ^ (O? 1) is a mixture of the two. 

We make the following modeling hypotheses. 

(Al) The displacements are small. Therefore, we state the problem in Lagrangian 
coordinates, in which the mass conservation is equivalent to the condition of a 
constant mass density Qq > 0. 

(A2) The substance is compressible, and the speed of sound does not depend on the 
phase. 

(A3) The evolution is slow, and we neglect shear viscosity and inertia effects. 
(A4) We neglect shear stresses. 

In agreement with (Al), we define the strain e as an element of the space T^y^ of 
symmetric tensors by the formula 

6 = V,u:=^(Vu+(Vu)^). (2.1) 
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Let 6 G TgyJI denote the Kronecker tensor. By (A4), the elasticity matrix A has the 
form 

Ae = X{e:5)S, (2.2) 

where " : " is the canonical scalar product in T^y^ , and A > is the Lame constant 
(or bulk elasticity modulus), which we assume to be independent of x by virtue of 
(A2). Note that A is related to the speed of sound vq by the formula vq = \JX/ Qq. 

We want to model the situation where the specific volume Vsoiid of the solid phase 
is larger than the specific volume Vuquid of the liquid phase. Considering the liquid 
phase as the reference state, we introduce the dimensionless phase expansion coefficient 
= {ysoUd — Viiquid)/Viiquid > 0, and we define the phase expansion strain e by 

-six) = ^{1 - X)S . (2.3) 

We fix positive constants Cq (specific heat), Lq (latent heat), 6^ (freezing point at 
standard atmospheric pressure), (3 (thermal expansion coefficient), and consider the 
specific free energy / in the form 

/ = Coe(l-\og(^))+^{{e-e{x)):Sy-l{e-e^)e:6 (2.4) 
V V6'c// 2go go 

where I is the indicator function of the interval [0, 1] . The stress tensor cr is decom- 
posed into the sum cr^ + <t^ of the viscous component cr*' and elastic component cr'^ . 
The state functions cr^' ,(t'^ , s (specific entropy), and e (specific internal energy) are 



given by the formulas 

a" = i^{et:S)S (2.5) 

(t' = go^ = i\ie:5-ail-x))-f3ie-e,))S, (2.6) 



e = f + es = Coe + —ie:S-ail-x)f + -0ee:S + Loix + nx)),i'2.8) 

2Qo Qo 

where u > is the volume viscosity coefficient. The scalar quantity 

p := -ust :S-X{e:S-a{l- x)) + [3(9 - 6,) (2.9) 

is the pressure, and the stress has the form cr = —p5. The process is governed by the 
balance equations 

— diver = f„o; (mechanical equilibrium) (2.10) 

fJoCf + divq = cr : £t (energy balance) (2-11) 
~loXt £ dxf (phase relaxation law) (2-12) 
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where 70 is the phase relaxation coefficient, is the partial sub differential with 
respect to Xi ^voi is a given volume force density (the gravity force) 



'■vol 



-QogSs, (2.13) 



with standard gravity g and vector <53 = (0,0, 1), and q is the heat flux vector that 
we assume in the form 

q = -kVO (2.14) 

with a constant heat conductivity n > 0. The equilibrium equation (12.101) can be 
rewritten in the form Vp = —Qog, hence 

p{x,t) = P{t)- gogx3 (2.15) 

with a function P of time only, which is to be determined. On dQ, we assume 
boundary conditions in the form 

u = (2.16) 

q-n = h{x){e-er), (2.17) 

with a given positive measurable function h (heat transfer coefficient), and a constant 
> (external temperature). Identity (12.161) means that the boundary is rigid. 
Other possibilities (elastic or elastoplastic boundary response) have been considered 
in another context ([ZllH]). 

By Gauss' Theorem, we have 

/ divu(x,t) dx = (2.18) 
Jn 

We have e : d = divu. Using (12. 9p . we write the mechanical equilibrium equation 
fl2J[5ll as 

udiv ut + A(divu - a{l - x)) - I3{0 - Be) + Pit) - g^g X3 = . (2.19) 
Integrating over Q and using (12.181) we obtain 

I I t/ I I t/ I I t/ 

We see that in liquid (x = 1) and at temperature 9 = 9c, the pressure p{x,t) vanishes 
on the "midsurface" of Q given by the equation x^lQl = J^x'^dx' . Hence, p{x,t) 
can be interpreted as the difference between the absolute pressure and the standard 
pressure. This difference is higher below and lower above the midsurface. 

Eq. (12.201) enables us to eliminate P(t) and rewrite (12.191) in the form 

z/divut + A(divu-a(l-x))-/3(^-^c) + 7^ / {aX{l - x) + Pi9 - 9^)) dx' 

l"l Jn 

= ^09(^3--^ I x'gdx') . (2.21) 
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As a consequence of fl2.4l) . the energy balance and the phase relaxation equation in 
fim]) - fl^:T^ have the form 



QqCqOi — i%/S.6 = ^{divutY — pOdivUf — [aX{diYu — a{l — x)) + QoLo)xt , (2.22) 

6 



QoloXt e aX{divu-a{l-x)) + QoLo{l--^ + dI{x)] , (2-23) 



■'c 



where d denotes the subdifferential. For simplicity, we now set 

c := QoCo , 7 := ^"070 , L := QqLq . (2.24) 

The system now completely decouples. For the unknown absolute temperature 9 , local 
relative volume increment U = divu, and liquid fraction x, we have the evolution 
system (note that mathematically, dl{x) is the same as Ldl{x) ) 

cOt-nM = uU^ - f3eUt- {aX{U -a{l-x)) + L)xt, (2.25) 



(2.26) 



1 



(aA(l -x)+ - Oc)) dx 



-IXt e a\{U-a{l-x)) + L(^l-pj+dI{x), (2.27) 

with boundary condition (12.171) . fl2.14p . that is, 

KVe ■n + h{x){e - Or) = 0. (2.28) 

We then find u as a solution to the equation divu = U in Q, u = on dQ. It 
is indeed not unique, and due to our hypotheses (A3), (A4), we lose any control on 
possible volume preserving turbulences. This, however, has no influence on the system 
(I2.25p - fi2.27p . which is the subject of our interest here. 

Let us describe the set of all possible stationary states. It follows from f l2.25p 
and (12.281) that the only temperature equilibrium is 9 = 9y ■ The equilibrium values 
Xoo and f/00 satisfy the system 



A(t/oo - a(l - Xoo)) = --pi y (1 - Xoo)dx' + ^ofi- (^3:3 - Xgdx' 



(2.29) 



-A([/oo-«(l-Xoo)) e ^ (^1 - + 9J(xoo) , (2.30) 
almost everywhere in fi, that is, 
^ (^-1) +P[£(l-^-)dx'-f?o^7(^X3- p^^dx^ g9/(xoo(x)). (2.31) 
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We claim that unlike in the case without gravity, (12.311) determines the equilibria 
uniquely. The set Q is connected. We can therefore define (a, 6) C M as the maximal 
interval such that Q fl (M^ x {0:3}) 7^ for X3 G (a, b) , and set 
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x'^dx' e {a,b) , i = b-a. (2.32) 



We introduce two dimensionless constants 



d = —— , G(] 



For water, we have for instance d ~ 0.055, Gq ~ 2 • lO'^/i. Eq. (12.311) then reads 

~ ^ M ~ + ^^"^ ~ "^'^ ^ ^^(Xoo(x)) . (2.33) 

The quantity 

is independent of x , so that the left hand side of (12.331) is positive for X3 < m + GoiZ 
and negative for X3 > m + GoiZ . By definition of the sub differential, we necessarily 
have 

f 1 if X3 < m + GoiZ , 
Xoo X = <^ ^ ^ 2.35 

[0 if X3 > m + GoiZ . 

Let Q{r) denote the set {x G : X3 > r} for r G M. Eq. (I2.35P states that the 
set f2(m + GoiZ) corresponds to the solid domain. We have |f^(r)| = for r > b, 
|n(r)| = for r < a, and 



1^ f {l-Xoo)dx' = F{Z) :- 

I I t/ 



\n{m + GoiZ)\ 



(2.36) 



Z = U^-1)+F(Z), (2.37) 



We easily identify Z as the only solution to the equation 

1 M 

d W 

since F is nonincreasing. We see that one of the following three cases necessarily 
occurs: 

(i) Z < {a — m) / {Goi) 1 F{Z) = 1. Then Xooi^) = a.e. in Q and we have pure 
solid with temperatures 

9r = 9,{1 + d{Z - 1)) <e,(^l-d(^l + 

(ii) Z > {b — m)/{Goi), F{Z) = 0. Then Xoo{x) = 1 a.e. in fl and we have pure 
liquid with temperatures 

Or = e,{i + dz) > eji + d-^ 
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(iii) (a - m)/(Go^) < Z < {b - m)/{Goi) , < F{Z) < 1 . Then Xoo{x) = a.e. in 
f2(m + Go£Z) , Xoo{x) = 1 a.e. in f2 \ Q{m + GqIZ) , and 

9,{l-d{l + — — ]]<9r<9jl + d 



Goi J J V GqI 

We observe that sohdification may take place at temperatures shghtly above 6c- For 
water in a container of £ = 50 cm height, the relative size of the "overheated ice tem- 
perature domain" is smaller than cI/Gq ^ 1.4-10"^, hence it is far beyond the standard 
measurement accuracy. The overheating is due to the fact the pressure decreases from 
the bottom to the top, as pointed out after formula f l2.20p . A quantitative charac- 
terization of this phenomenon is given by the so-called Clausius-Clapeyron equation, 
which relates the freezing temperature with the pressure. It can be derived here as 
follows. The equilibrium relative pressure p^o depends only on x^, and is given, by 
virtue of fl2J[5D and (E^Ql), by the formula 

Poo(x3) = ^ / (1 - ^°o) dx' + (3{er -dc) + ^J x', dx' - Qogxs . (2.38) 

I I t/ I I •/ 



Using (I2.3ip . we obtain 
L fOr 



P{Or - 9c) + Poo(x3) e dI{xoo{x)) . (2.39) 



The phase interface at temperature is located at level 0:3 if the right hand side of 
fl2.39p vanishes. Setting Lp = Lq — a[36c/ Qo , we thus obtain the Clausius-Clapeyron 
condition for phase transition in the form of [IDl Equation (288)], namely 

— 9c a^c (Viiguid — Vsoiid)9c 
In terms of the new variables 6,U,x, the energy e and entropy s can be writ- 



ten as 



coO + -^{U - ail - x) f + -OcU + Lo(x + /(x)) , (2.41) 
cologf^)+^X+^t/, (2.42) 



and the energy and entropy balance equations now read 
d 



dt 



Qoe{x,t) — QogxsU{x,t)) dx = / h{x){6r — 6) ds{x) , (2.43) 



an 



Qost + div^ = -L^ + lx'^ + -Uf > 0, (2.44) 



A f g^s{x,t)dx = I ^{er-e)ds{x) (2.45) 



, ,'fi:|V6'P 7 2 ^rrz \ 1 
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The entropy balance (I2.44p says that the entropy production on the right hand side is 
nonnegative in agreement with the second principle of thermodynamics. The system 
is not closed, and the energy supply through the boundary is given by the right hand 
side of f lCTD . 

We prescribe the initial conditions 

e{x,o) = e\x) (2.46) 

f/(x,0) = f/°(x) (2.47) 
X{x,0) = x\x) (2.48) 

for a; G fi, such that f^Uo{x) dx = 0, and compute from f l2.4ip - fl2.42p the corre- 
sponding initial values e° and for specific energy and entropy, respectively. Let 
= f^QQe'^dx, 5"° = J^Qos'^dx denote the total initial energy and entropy, re- 
spectively. From the energy end entropy balance equations fl2.43p . fl2.45l) and using 
the condition J^Udx = 0, we derive the following crucial (formal for the moment) 
balance equation for the "extended" energy go{e — 9rs): 

c9 + ^{U - a{l - x)? + Lx- QogxsU^ {x,t) dx (2.49) 



+ / / ^{9r-enx,T)ds{x)dT 
Jan ti 



= E'- ^r5° - Qog ^ XsU'ix) + ^ (clog (^^j + ^x) {x, t) dx . 

We have log(^/^c) = log(^/2^r) -log(^c/20r) < (^/2^r) - 1 -log(^e/2^r) , hence there 
exists a constant C > independent of t such that for all t > we have 

Jje + u') (x, t) dx + J^J^ + y + t) ^^-^^^ 
^{er-e)\x,T)ds{x)dT < c. 



Jdn 



3 Existence and uniqueness of solutions 

We construct the solution of f l2.26p - fl2.27p by the Banach contraction argument. The 
method of proof is independent of the actual values of the material constants, and we 
choose for simplicity 

L = 2, g = c = qq = 9c = a = (3 = = K = \ = V = \ . (3.1) 



8 



System fl2.25l) - fl2.27l) with boundary condition fl2.17p then reads 
I dtw{x)dx+ I Ve-Vw{x)dx = I (U^ -eUt- (u + x + l)x?)w{x)dx (3.2) 



h{x){e - er)w{x) ds{x) 
an 



Ut + U = -x + »+(i'3-|^/^4di-'| (3.3) 



Xt + U + x + dI{x) 3 26 - 1, (3.4) 

where (13.21) is to be satisfied for all test functions w G W^''^{Q) and a.e. t > 0, while 
(I3.3l) - (l3.4p are supposed to hold a.e. in Qoo := f2 x (0, oo) . 

In this section we prove the following existence and uniqueness result. 

Theorem 3.1 Let < 6^^ < 6^ < 6* and pq ^M. be given constants, and let the data 
satisfy the conditions 

00 e W^'^{Q)r\L°°{n) , 6, < e^{x) < e* a.e., 

e L'^^iSl), /^?7°(x)dx = 0, < x°(a;) < 1 a.e. 

Then there exists a unique solution {6, U, x) to ^K2)-(3^, \2.4^6\} - [2.4-9\) , such that 
e > a.e., X ^ [0,1] a.e., U,Ut,Xu 0,1/9 E i:°°(l^oo), Ot, A9 E L'^{Vt^), and VO E 
L-(0,T;L2(fi))nL2(fi^). 

Remark 3.2 For existence and uniqueness alone, we might allow the external tem- 
perature Or to depend on x and t, and assume only that it belongs to the space 
Wl^l{Q, oo; L'^{dVt)) n L-^^iOVt x (0, oo)). For the global bounds, the assumption that 
be constant plays a substantial role. 

The proof of Theorem 13.11 will be carried out in the following subsections. Notice 
first that the term — 9Ut — {U + X + ^)Xt on the right hand side of (13. 2p can be 
rewritten alternatively, using (13. 4p and (13. 3p . as 



ut-eut-{u + x + i)xt = ut -eut + xi-'^oxt (3.5) 



ix + U)Ut+[x^-^^ l^x',dx']Ut 
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{-X + 0) dx')Ut-iU + x + l)Xt 
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We now fix some constant R > and construct the solution for tlie truncated system 

[ etwix)dx+ I VB-Vw{x)dx = I {U'{ + xl-QR{d){Ut + 2xt))w{x)dx{'i.Q) 
Ja Jn Jn 

-[ h{x){9-9r)w{x)ds{x) WweW^^^Q), 
Jan 

Ut + U = -X + QBiO)+(^xs-^J^x',dx'^ (3.7) 

I I t/ 

Xt + U + x + dI{x) ^ 2Qr{9)-1 (3.8) 

first in a bounded domain Qt := x (0,T) for any given T > 0, wliere Qr is tlie 
cutoff function Qr{z) = min{z~^,R}. We tlien derive upper and lower bounds for 
6 independent of R and T, so that the local solution of fl3.6l) - fl3.8p is also a global 
solution of fl3.2l) - fl3.4p if R is sufficiently large. 



3.1 A gradient flow 

System fl2.26p -( !2.27p can be considered as a gradient flow similarly as in the case 
without gravity and with an elastic boundary, see [71 Section 4.1]. Set 



U 
X 



(3.9) 



Hv) = ^ (^(^ - (1 - X)f + 2x (1 - ^r) - xsU + /(x)) dx (3.10) 
+ / f^dx {l-x + X3)dx + C^, 



with a constant > to ensure that ipi'^) ^ 0. Every solution U of (12.260 neces- 
sarily satisfies the condition J^U dx = , hence fl2.27p can be written in the form 

Xt + U + x+l-2e--^ I f/dx + <9/(x) 9 0, (3.12) 

I I £^ 

and (12.260 . fl3.12p are in turn equivalent to the gradient flow 

v + d4){v)3f (3.13) 

in L^(0) X L'^{Vt) . We have used the obvious identity 9 — 1 — J^{9 — 1) dx = 
{9 — 9r) — 1 / \Q\ J^{9 — 9r) dx . We state here the following Lemma, whose proof can 
be found in [3 Lemma 4.3]. 
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Lemma 3.3 Let f,f belong to L^(0, oo;i/). Then lim(_^oo 'i^(^) = 0. 



We apply the above result to the case H = L^(f2) x L^(f2) , and v , f , ip as above 
and we see that Eqs. fl3.7l) -( !3l8l) with 6 replaced by 9 can be equivalently written as 
a gradient flow fl3.9p - (l3.13l) . For its solutions, we prove the following result. 

Proposition 3.4 Let the hypotheses of Theorem \3.1\ hold, and let a function 9 G 
Ll^^{{),oo;L'^{VL)) he given. Let (t/, x) he the solution of [3^} - [3T^} . Then there 
exists a constant Cq, independent of x,t and R, such that a.e. in fioo have 

\U{x,t)\ + \Ut{x,t)\ + \xt{x,t)\ < Co(l + . (3.14) 

Let furthermore 9i,92 G Lf^^{0, oo; L'^{Q)) be two functions, and let (f/i, Xi); (^2, X2) 
be the corresponding solutions of 1^3. 9\) - [3713) . Then the differences 9d = ^1 — ^2? 
Ud = Ui — U2 , Xd = Xi ~ X2 satisfy for every t >0 and a.e. x G the inequality 

j^md)t\ + \{Xd)t\){x,T)dT<C,{l+t) (\Ux,T)\+t\9d{T)\2) dr, (3.15) 

where the symbol | ■ I2 stands for the norm in L'^{Q) . 



In what follows, we denote by Ci,C2,... any constant independent of x,t 
and R. 

Proof. We rewrite f l3.9p -( l3rT3l) as two scalar gradient flows 

Ut + dtpi{U) = a, (3.16) 
Xt + dij2{x) ^ b, (3.17) 

where MU) = \U\ 1P2 = \x^ + /(x), a = Qr{9) - x - ^L(Qi?(^) - x) dx + 
(•^3 ~ ^ /o^sd^^'j , h = 2Qfi{9) — 1 — U . The bounds (13.141) are obvious. To prove 
fl3.15p . we consider two different inputs. As above, we denote the differences {}i — {}2 
by {}d for all symbols {}. By [HI Theorem 1.12], we have for all t > and a.e. x G f2 
that ^ ^ 

' md)t\ + \{Xd)t\){x,T)dT<2 [ {\ad\ + \bd\){x,T)dT. (3.18) 

Jo 

We multiply the difference of (I3.16p by Ud, the difference of (13.170 by Xd, and sum 
them up to obtain that 

(f/d)it/d + (Xd)ad + (?7rf + Xd)'< |^d|(|f/d| + 2|xd|)-(^p^(g«(^)d-Xrf)dx^ Ud a.e. 

(3.19) 

We first integrate fl3.19p over Q. Using the symbol | ■ I2 for the norm in L'^{Q), we 
get for a.e. t > that 

{\Ud\l + \Xd\l) < \9M\Ud\2 + 2\xd\2) < V5\9d\2 {\Ud\l + IXdllY^'' ■ (3.20) 
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Hence, -§f.{\Ud\2 + IXdli)"*^^^ — v^l^dh a.e., and integrating over t, we find that 

{ml + \Xd\lY^\t) < V5 f\9,{T)\,dT. (3.21) 
Using again fl3.19p . we find for a.e. (x, t) G Qoo the inequahty 
l^t^U! + xl) {x,t)<Mx,t)m+2\xMx,t) 

+ Ci|f/,|(x,t) (Mt)\2+ t \e,{x,T)\2dT) . (3.22) 



Hence, we get 

Wt + ^"^'^^ - + Xlf'\x,t) {\e,\{x,t) + |^,(t)b + \Ur)\2dr 

(3.23) 

which in turn imphes that 

{\Ud? + \Xd?y'^ (x, t) < C3 (^^* |^,(x, r)| dr + (1 + t) |^,(r)|2 dr^ a.e. (3.24) 
This enables us to estimate the right hand side of fl3.18p and obtain the bound 

md)t\ + \{Xd)t\){x,T)dT<Ci [ ({l+t)\ed\{x,T)+t{l + t)\UT)\2) dr (3.25) 



'0 ^0 ^ ' 

for a.e. x and all t>0. This completes the proof. ■ 

3.2 Existence of solutions for the truncated problem 

We construct the solution of fl3.6p -( !3l8l) for every i? > by the Banach contraction 
argument on a fixed time interval (0, T) . 

Lemma 3.5 Let the hypotheses of Theorem \3.1\ hold, and let T > and R > 
be given. Then there exists a unique solution {9, U, x) to li3.6\) - (37^) . ^2.46\ )- [2.48^ , 



such that U G W^'^inr) , > a.e., Xt, 0,1/0 G L^iflr), 0t,A9 G L'^iflr), and 
VO G L°°(0,r; L2(fi)). 

Proof. Let 9 G L^(f2r) be a given function, and consider the system 

/ Otw{x)dx+ [ \/0-\/wix)dx = [ (Uf + Xt -QRiO)iUt + 2xt))wix)dx (3.26) 
Jn Jn Jn ^ ' 

~[ h{x){0-Or)w{x)ds{x) ^weW^^^n), 
Jdn 

Ut + U = -x + Qn{0)+(^xs-^J^x',dx''^ (3.27) 

-^Jj-x + QRiO))dx', 

Xt + U + x + dI{x) 3 2QR{e)-l (3.28) 
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Equations fl3.27l) - fl3.28l) are solved as a gradient flow problem from Subsection 13.11 
while fl3.26l) is a simple linear parabolic equation for 6. Testing fl3.26p by ^j, we 
obtain by Proposition 13.41 that 

I [ e^dxdt + supess ( [ \\/9\^dx+ [ h{x){9 - 9rf ds{x)] (3.29) 
Jo Jn ie{o,T) \Jn Jan J 

< T\n\ (Co(l + i?)(2Co(l + R) + 3R)f =: Mr . 

Hence, we can define the mapping that with 6 associates the solution 6 of fl3.26l) - 
(I3.28P with initial conditions fl2.46l) - fl2.48l) . We now show that it is a contraction on 
the set 

Et,r := {9 e L^inr) : conditions fICTD - flOill hold} , (3.30) 

where 



V9 G L°°(0,T; L\n))- 



(3.31) 
(3.32) 



T 



Jn 



te{o,T) \Jn 



^^dxdt + supess / |V^pdx+ / h{x){9 - 9rY ds{x) ) < Mr; (3.33) 



an 



9{x,0) = 9%x) a.e. 



(3.34) 



Let 9i,92 be two functions in St,_r, and let {9i,Ui,xi) , (6*2, f^2, X2) , be the corre- 
sponding solutions to fl3.26p - fl3.28p with the same initial conditions 9^,U^,x'^- We 
see from fl3.29|] that 9i,92 belong to 'Et,r- Integrating Eq. fl3.26p for 9i and 6*2 with 
respect to time and testing their difference by w = 9^ := 9i — 92, v/e obtain, using 
Proposition 13.41 that 



9l{x,t) dx + — 



9d{x,T) dr 



da; 



h{x) 



an 



9Jx, t) dr 



ds{x) 



md)t\ + \iXd)t\ + mix,T)dT)9dix,t)dx a.e. (3.35) 



From fl3.15p and Minkowski's inequality, it follows that 



md)t\ + \{xMi;r)dr 



< Ce{l + tf / |^,(r)|2dr 



< Ce{l + tf{t / |^d(r)|2dr 



1/2 



By Young's inequality, we rewrite (13.351) as 



V / 9d{x,T)dr 



9l{x,t) dx + — 



< Cr{l + R'){l + tf |^d(r)|^dr a.e. 







\ 


dx + h{x) 


/ 9d{x, r) dr 


ds{x) j 


Jan 


Jo 




a.e. 




(3.36) 
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Set e2(t) = /o*|^d(r)|2dr, e^{t) = /q |^d(r)|i dr . Integrating (Km with respect to 
time, we obtain 

e''{t) <Cr{l + R^) [ (l + r)^e2(r)dr. (3.37) 
Jo 

We set Cr := {Cjil + R^)/Q) and introduce in L°°(0,T) the norm 

||w;||c := sup e'^^^^+^^VWI • 
Te[o,T] 

Then ||e||^ < |||e||^, and hence the mapping ^ i— >■ ^ is a contraction in L'^{Qt) with 
respect to the norm induced by || ■ ||c. The set St,_r is a closed subset of L'^{Qt)- 
This imphes the existence of a fixed point 6 G H^ /j, which is indeed a solution to 
(I3.6l) - fl3.8p . The positive upper and lower bounds for 6 follow from the maximum 
principle (for the proof cf. [TJ Section 4.2]). This completes the proof of Lemma [3.5[ 



3.3 Proof of Theorem [3lT] 

The unique solution {6, U, x) to (l3J])-(l3Sj), ^M)-^M) exists globally in the whole 
domain Qoo- We now derive uniform bounds independent of t and R. Take first for 
instance any R > 26* . We know that the solution component 6 of f l3.6p -( l3l8l) remains 
smaller than i? in a nondegenerate interval (0,T) with T > e*/{CQ{l + Rf). Let 
(0,To) be the maximal interval in which 6 is bounded by R. Then, in (0,To), the 
solution given by Lemma 13.51 is also a solution of the original problem (I3.2l) - fl3.4l) . 
Moreover, due to estimate (12.501) . we know that 6 admits a bound in L°°{0, Tq; L^{Q)) 
independent oi R. In order to prove that Tq = +oo if R is sufficiently large, we need 
the following variant of the Moser iteration lemma, whose proof can be found in [71 
Prop. 4.6]. 

Proposition 3.6 Let Q C be a bounded domain with Lipschitzian boundary. 
Given nonnegative functions h G L^{dQ) and r G L°°(0, cxo; L'^(f2)) with a fixed 
q > N/2, |r|ioo(o,oo;L9(n)) =: r* , an initial condition G L°°{Q) , and a boundary 
datum vr G L°°{dfl x (0, oo)), consider the problem 

Vt-Av + v = r{x,t)n[v] a.e. m n X (0, cx)) , (3.38) 
Vf ■ n = —h{x){f{x,t,v{x,t)) — vr{x,t)) a. e. on dO, x {0, oo) , (3.39) 
v{x,0) = v° a.e. in Q , (3.40) 

under the assumption that there exist positive constants m, Hq, Cf,V,Vr,EQ such that 
the following holds: 

(i) The mapping Ti : L'^^iVL x (0,oo)) L'^^iVt x (0,oo)) satisfies for every v G 
L'^ci^ ^ (0, cxd)) and a.e. {x,t) eQx (0, oo) the inequality 

v{x,t)n[v]{x,t) < Ho\v{x,t)\ (^1 + \v{x,t) \ + i{t-T)\v{x,T)\dT 
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where ^ G H^^'"'^(0, oo) is a given nonnegative function such that 

m<-mm (3.41) 

(ii) f is a Caratheodory function on Q x (0, oo) x M such that f{x,t,v)v > CfV^ 
a.e. for a// f G M. 

(iii) |i;°(x)| <V a.e. in Q . 

(iv) \vr{x,t)\ < Vp a.e. on dil x (0, oo) . 

(v) System 38\ )- ^y.40\ ) admits a solution 

V G <f (0, oo; {W'^m n LlM oo; w''\n)) n LZi^ X (0, oo)) 
satisfying the estimate 

\v{x,t)\dx < Eq a.e. in {0, oo) . 

Then there exists a positive constant C* depending only on |/i|Li(en) > C/? Hq such 
that 

\v{t)\L^(^) < C*m&x{l,V,Vr,Eo} for a.e. t > 0. (3.42) 

We now finish the proof of Theorem 13.11 by showing that To introduced at the 
beginning of this subsection is +oo if R is sufficiently large. Using 03.31) . we obtain 
that 

\U{x,t)\ < Cs(l + e^-'e{x,T)dT^ a.e., (3.43) 

\Ut{x,t)\ < Cg(^l + e{x,t) + e^-'e{x,T)dT^ a.e., (3.44) 
hence also (cf. (13. 4p ) 



\Xt{x,t)\ < Cw\^l + e{x,t) + J e^~'9{x,T)dTj a.e. (3.45) 

By (12.501) . the function U is in L°°(0, oo; T^(f^)) and the bound does not depend on R. 
Eq. (13. 2p . with 6 added to both the left and the right hand side, thus satisfies the 
hypotheses of Proposition 13.61 for = 3 and q = 2. This enables us to conclude that 
6{x,t) is uniformly bounded from above by a constant, independently of i?, so that 
9 never reaches the value i? if i? is sufficiently large, which we wanted to prove. By 
(I3.43p - (l3.45p . also U , Ut, and Xt are uniformly bounded by a constant. 

We proceed similarly to prove a uniform positive lower bound for 6 . Set Rq := 
sup^, and in Eq. (13. 6p with R > Rq put w = —w/9, w G W^''^{fl). For a new 
(nonnegative) variable v{x,t) := logi?o — log6'(x, t) we obtain the equation 

Vtw{x)dx+ I Vf ■ Vty(x) dx + I h{x) [ — i\ w{x) ds{x) (3.46) 
Jn Jan 



9 ^2 
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9 

+ Ut + 2xt ] w{x) dx . 



We now set 

n[v] = sign{v) [-^^^ - ^ + ^* + 

and check that the hypotheses of Proposition [H^ are satisfied with f{v) = {9r/Ro){G^ — 
1) , vr = (Ro — Or)/Ro, r = 1, and vT-C[v] < 2C\v\, where C is a common upper bound 
for Ut and Xt ■ Hence, v is bounded above by some v* , which entails 9 > Roe~^* . 
This concludes the proof of Theorem I3.1[ 



4 Long time behavior 

We have the following statement. 



Proposition 4.1 Let the hypotheses of Theorem \3. 1\ hold. Then we have 

[ ( [ {el + f/2 + xl + |V^n dx + / h{x){e - Ovf ds{x)] dt < oo, (4.1) 
Jo \Jn Jdn J 

\im ( [ (U^ + Xt + \^On{x,t)dx+ [ h{x){e -ery{x,t)ds{x)] = 0. (4.2) 
\Jn Jan J 

Furthermore, letting t tend to oo , the temperature 9 converges strongly in W^''^{VL) to 
its equilibrium value 9^ , and also both x(x, t) and U{x,t) converge strongly in L^{Q) 
(hence, strongly in every Lp{Q) for p < oo) to their respective unique equilibrium 
values Xoo and Uoo defined in Section\^ 



Proof. The proof of the relations fl4.ip - fl4.2l) exactly follows the argument of [71 Propo- 
sition 5.1]. For the convergence of the whole trajectory, it seems easier to prove it 
directly without referring to the general theory of dynamical systems with a unique 
equilibrium (see, e.g., [21 Lemma 2.4]). 

We eliminate from f l2.26p -( l2.27l) the term U — a{l — x) ^"^^ obtain, using the 
notation (12321) - dOil) . the inclusion 

A{x, t) + p ^(xoo - X) dx' + Qog{m + Go^Z - xg) G dl{x) , (4-3) 
where we set 

A{x,t) ■.= uU,-\t+(^-/3)i9-9r) + ^ [ i9-9r)dx'. 
a \a9c J \VL\ J^ 

From (14.31) it follows that 
T^(Xoo -X) I (Xoo - x) dx' + ^'ofi'(Xoo - x)("^ + Go^Z - xs) < (xoo - X)A{^. t) (4.4) 

II t/ 
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for a.e. (x, t) G floo ■ Integrating fl4.4p over fl and using (12.351) we obtain 



W\ 



/ (Xoo-x)da;' + Qog / \xoc>-x\\m + Go£Z -x'^\dx' < / |xoo-xll^(a^',i^)l dx' . 
Jn Jn Jn 

(4.5) 

By virtue of (14.21) . we have limt^oo /q |A(a;', t)| dx' = 0. Hence, the right hand side 
of (14. 5 p converges to as t ^ oo. Since x is a priori bounded, the weighted 
convergence of Xoo — X with weight |m + GqEZ — x-sl imphes that x ~^ Xoo strongly 
in L^{Q) as t ^ oo. The convergence U Uoo follows directly from (I2.26p . (12.291) . 
and m 
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